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We discuss the high temperature superconductors in a regime where the antiferromagnetic (AF) 
correlation length is only a couple of lattice spacings. In the model proposed here, these short-range 
AF fluctuations play an essential role in the dressing of the carriers, but the attraction needed for 
superconductivity (SC) arises from a transverse phonon oxygen mode with a finite buckling angle 
as it appears in YBa2Cu307-i. A simple fermion-phonon model analog to the Holstein model is 
introduced to account for this effect. We argue that the model has a d^2_y2-wave superconducting 
groundstate. The critical temperature (Tc) and the O-isotope effect coefficient (ao) vs hole density 
(x) are in qualitative agreement with experiments for the cuprates. The minimum (maximum) of 
cto (Tc) at optimal doping is caused by a large peak in the density of states of holes dressed by AF 
fluctuations, as discussed in previous van Hove scenarios. 
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Since the discovery of high temperature superconduc- 
tors, the origin of their pairing mechanism has been con- 
troversial. Numerous studies have shown that normal 
state properties deviate from a conventional Fermi liquid, 
and as a possible explanation several authors proposed 
the AF correlations as responsible for such nonstandard 
behavior. Theories based on AF pairing mechanisms an- 
alyzed using diagrammatic and numerical techniques pre- 
dict da;2_j,2 SC. j^] Josephson junction experiments and 
angle-resolved photoemission (ARPES) data are consis- 
tent with such d^2_y2 condensate. [|^,^ 

However, there are still some problems with this ap- 
proach. For example, the AF correlation length, £,af, in 
the normal state of the high-Tc cuprates at optimal dop- 
ing, Xopt, niay not be robust enough to induce SC. NMR 
measurements in YBa2Cu307_5 suggest ^af/cl ~ 2 — 3 
Q (a is the Cu-Cu lattice spacing). Inelastic neutron 
scattering studies for the same compound only show a 
broad peak at q = (vr, tt) in the dynamical spin struc- 
ture factor, g Weak "shadow bands" in Bi2212 ARPES 
data at room temperature have been interpreted as pro- 
duced by short-range AF correlations (although superlat- 
tice effects have not been ruled out), Numerical stud- 
ies have shown that these weak AF-induced bands are 
quantitatively reproduced by the doped t-J model in a 
regime where ^af/o, ~ 2 — 3, in agreement with NMR. 
While many numerical studies suggest that this appar- 
ently small AF-correlation can nevertheless substantially 
affect the quasiparticle (q.p.) dispersion (one-particle 
Green's function), ||] it has not been shown that it can 
also lead to pairing (two-particle Green's function, with 
hole Coulombic repulsion included) as it may occur in 
the regime ^af/o- > 1. 

A second problem for electronic mechanisms in general 
is the O-isotope effect observed in the cuprates. While 
experiments have shown that the coefficient ao is very 
small at Xopt, it increases in the underdoped and over- 
doped regimes reaching values comparable to the BCS 



limit aecs — 0.5. |S[| It may be argued that this effect 
is caused by spurious changes in x after the replacement 
O^^ O^^, but until a realistic calculation proves it, the 
experimental data cannot be simply neglected. 

In this paper we assume that the normal state £,af 
at Xopt is not large enough to produce magnon mediated 
pairing, and thus we discuss possible alternative ideas 
that may explain SC in the cuprates. We argue that a 
small ^AF can still be important for the normal state hole 
dispersion since carriers are much affected by the sur- 
rounding local spin environment. Short AF correlations 
can modify the hole dispersion reducing the bandwidth 
and producing anomalous "flat-bands" as observed in the 
t-J and Hubbard models, ||| and in ARPES data, g The 
flat-bands induce a robust peak in the density of states 
(DOS) at the top of the valence band that boost Tc and 
produce an "optimal doping" when the chemical poten- 
tial reaches the peak, once a source of hole attraction 
exists. This combination of the AF and van Hove sce- 
narios pO[ | is similar in spirit to previously discussed van 
Hove theories. The key difference is the origin of the 
DOS large peak which in Ref. ||l^ was attributed to AF 
correlations. In the present paper, pairing is caused by 
a p/ionoTiic-induced attraction supplemented by a hole 
dispersion modulated by short distance AF correlations 
that tend to prevent double occupancy and favor intra- 
sublattice hopping. Thus concepts of both phononic and 
electronic theories are here mixed in a single scenario. 

However, an immediate problem with this idea is the 
symmetry of the SC condensate. Evidence is accumu- 
lating in favor of d,j.2_y2 SC which is natural in elec- 
tronic AF-theories, but seems unnatural in conventional 
phononic theories. For example, the Holstein model cou- 
ples electrons to on-site oscillators ||T^ leading to a uni- 
form s-wave condensate. To obtain d-wave phononic 
SC, we introduce a modification of the Holstein model 
with oscillators located at the oxygens of a two dimen- 
sional (2D) square lattice rather than at the copper 
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sites (Fig. la, b). While this modification seems "ad-hoc", 
physical realizations involving the buckling mode of oxy- 
gen in the cuprate exist, as shown below. We further ar- 
gue that the effective hole-hole interaction produced by 
this model favors d-wave SC at low carrier concentration 
once the AF induced hole dispersion is used. 
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FIG. 1. (a) Graphical representation of the standard Hoi- 
stein model with an oscillator attached to each lattice site i; 
(b) In the new model Eq.(l) the oscillators are attached to 
the lattice Cu-Cu links (oxygen positions); (c) Oxygen buck- 
ling mode studied in this paper, zq is the oxygen equilibrium 
position, and Zi^~, is the displacement from equilibrium in the 
7 direction. 

The modified Holstein model proposed here is 
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where eApO^) is the hole dispersion obtained after the 
influence of AF correlations has been taken into account 
(this is important since the resulting q.p.'s are weakly 
interacting and thus the peak in the DOS is not removed 
by correlations or disorder jl^). c{c') is a fermion de- 
struction (creation) operator located at Cu with standard 
fcrmionic anticommutation relations (these fermions rep- 
resent "holes" in the cuprates), •hi^ — c\^Cia is the num- 
ber operator, i labels sites of a 2D square lattice, 7 = x, y 
are unit vectors along the axis-directions, and Zi^-y and 
Pi^~i are the coordinate and momentum of the oscillators 
of mass M and frequency lo at half-distance between ad- 
jacent sites to mimic oxygens. As hole dispersion we 
use eApC^) / gV = Q.IG^ coskx cosky + Q.OAib {cos2kx + 
cos2ky), which corresponds to the hole dispersion in an 
AF background at half-filling, |^ but it should be ap- 
proximately valid also at finite x as long as ^af is not 
negligible. Holes move within the same sublattice. The 
spin index in Eq.(l) does not play an important role and 



working with spinless fermions leads to similar results. 
Using Eq.(l) we will study small fermionic (hole) densi- 
ties to mimic the physics of the cuprates. 

What cuprate phonons can lead to Eq.(l)? The in- 
plane breathing mode where O oscillates along the Cu- 
Cu link is not useful since it produces an effective hole 
n.n. repulsion. In principle 0-phonons transverse to the 
Cu-Cu link have a coupling quadratic in the displace- 
ment that is negligible. However, several authors ^3 14 
noticed that the buckling of the Cu-O-Cu link leads to 
a linear electron-phonon coupling in YBa2Cu307_5. In 
this tilting mode (Fig.lc) the 0-atom oscillates in the z- 
direction about an equilibrium position zq with a buck- 
ling angle /? = 5° and a frequency that here we take as 
fi^buck = 12 meV. 13| While buckling effects are not 
present in all cuprates, nevertheless it serves our pur- 
pose of identifying at least one phononic mode that in 
combination with strong correlations can lead to d-wave 
SC. Note that in this exploratory study we are also ne- 
glecting the coupling to other phonons modes of similar 
energy that may lead to repulsive interactions. 

Let us discuss recent work related to the ideas de- 
scribed here. Song and Annett ||l^ studied d-wave 
phononic SC using the in-plane 0-breathing mode. How- 
ever, further analysis showed that this mode does not 
lead to d-wave SC. ||l^ They also studied the 0-buckling 
mode with a tight binding dispersion finding d-wave SC, 
but concluded that it would not produce a large enough 
Tc- jl^ The key difference with our approach is that we 
here use an AF-induced DOS with a large peak which 
boots Tc to high values. In other related work, Yone- 



mitsu et al. |17| coupled the apical phonon modes to the 
t-J model forming a SC polaron pair condensate. AF cor- 
relations affect the hole propagation as in our approach, 
but they found "nodeless" d-wave or p-wave pairing, con- 
trary to our d,j.2_y2 result. Finally, using a mean-field ap- 
proximation for the t-J model with phonons, Normand et 
al. |lj] studied anomalies at Tc associated to the YBCO- 
buckling mode, ao was reported in agreement with ex- 
periment but the calculation was done only at Xopt ■ 

Returning to the main idea, let us derive the buckling- 
mode induced holc-phonon coupling. Consider the 
Coulomb energy of carriers at the Cu-ions in the pres- 
ence of the n.n. 0-ions, ||l8| 

_ ee* _t _ . 1 1 X 

^Coulomb - — 2^ "^-^'^^IRi-riJ + |Ri-ri _J^' 

(2) 

where e is the electron charge, e* = — 2e is the O- 
ion charge, e is the dielectric constant, Ri denotes 
the Cu-positions which are assumed non-fluctuating, 
and Ti -y denotes the vibrating 0-positions (interactions 
at distances larger than a are assumed negligible due 
to screening effects). In Eq.(2) we define ri__^ — 
Ti-rj^j. The Cu-0 distance can be expanded in the 
small 0-ion displacement in the z-direction, zi^-y, as. 
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I Ri - r;,^ 1= 1(1 + ^) + ^Zi^^ + ...,. The relation 
a ^ S> zi^ry is assumed. The hole-phonon interaction 
becomes 
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X]cLci<T(^i,7 + 2i,_^). (3) 



In k-space we arrive to the fermion-boson Hamiltonian, 

Hh-ph = ^J2 5q,7cLck-q<T(&q,7 + fc-q,^), (4) 
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where N is the number of Cu-sites, &q,7(&q-y) is the de- 
struction (creation) phonon operator with momentum q, 
and the electron-phonon couphng has the form 
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where M is the the 0-mass. ||T^ The strength of the cou- 
phng is estimated as ga=o.7 ~ (1 ~ 2) • 10~^ eV where 
we used e « 10-20 a w 3.8l, zq « 0.171 [fi, 
and M — 16a.u.. This is in agreement with Ref. [13|, 
and it is at least one order of magnitude less than the 
typical electron-phonon coupling strength in normal met- 
als, which is natural due to the reduction caused by the 
geometric factor zq/u. 

Although a detailed study of Eq. (4) would require the 
Eliashberg equations, here we only analyze the nonre- 
tarded version of the effective hole- hole phonon- mediated 
interaction. This simplification should not change the 
symmetry of the SC condensate and other qualitative 
features discussed below. Standard manipulations lead 
to the interaction (natural units): 
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where V^p^ is given by 
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Since we are considering a dilute gas of holes and a short- 
range potential Eq. (7) , this problem can be studied with 
the gap equation: Ap = - X;k(Kp^^k/2£;k) tanh ^, 

where = ^JleAF^ - m]^ + ^^^^ Ak is the SC 
gap. The numerical solution of this equation produces SC 
in the d^2_y2 channel. To obtain Tc vs x the linearized 
gap equation was solved. ||2l[] Naively, the small hole- 
phonon coupling induced by the buckling mode should 
produce a small T;.. However, the large peak in the hole 
DOS caused by fiat bands can boost T^. (lO[|l|] The re- 
sults are shown in Figs. 2 and 3. Tc is maximized when fi 
reaches the maximum in the DOS producing an "optimal 
doping". |ri| For the couplings and the approximations 



described in this paper, which are standard, Tc reaches 
30K which can be made higher by tuning the q.p. disper- 
sion to increase the boosting DOS. Thus, the experimen- 
tal symmetry of the SC phase, and the x-dependence of 
Tc can be qualitatively reproduced by the small fermion- 
phonon coupling induced by the buckling mode of YBCO, 
if e/ii?(k) is used as dispersion. 



40 



35 



TJK) 



30 



25 



20 



15 



10 




5 10 15 20 25 30 35 40 
X (%) 

FIG. 2. Tc vs X using e(k)AF as dispersion. The coupling 
gq=o,-y is set to 0.018 eV. The SC state is d-wave. 

Finding d-wave SC from Eqs.(6,7) is natural. To visu- 
alize this effect, assume interacting on-shell holes. Then, 
transforming Eq.(6) into real space we get 
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where (ij) denote n.n. sites. The first term in Eq.(8) 
is an on-site attraction which is suppresed trivially by 
the hard-core properties of the holes in the t-J model 
from which the dispersion is derived. |^ The second 
term provides the n.n. attraction that leads to an in- 
teraction of the form corresponding to the well-known 
"t-U-V" (U > 0, V < 0) model that has the tendency to 
form d-wave condensates. ||2^ However, there are sev- 
eral crucial differences between our model and the t-U-V 
model. More remarkable is the fact that the AF disper- 
sion used here allows the formation of a d-wave conden- 
sate at low particle density, while the t-U-V model (where 
a coskx -\- cosky dispersion is used) has s-wave SC in this 
regime [p2| . Optical conductivity measurements in un- 
derdoped cuprates clearly show that the number of car- 
riers grows like the number of holes, and thus we should 
study a dilute gas of quasiparticles as carried out in this 
paper. Our ideas go beyond previous studies of d-wave 
SC that have used the t-U-V model in spite of its short- 
comings like having a Tc maximum at half-filling and a 
strong competition with phase separation, p^ ] Combin- 
ing the potential Eq.(7) with the AF dispersion Q is the 
proper way to mimic the phenomenology of the cuprates. 
In addition, the phonons in Eq.(l) produce an isotope ef- 
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feet which does not exist in the electronic t-U-V model 
used before in the literature. Note also that the instan- 
taneous Coulombic repulsion at distance a would tend to 
suppress an effective attraction of electronic origin, like 
magnon mediated interactions, since it does not lead to 
a substantial retardation effect. On the other hand, the 
retardation intrinsic to phononic mechanisms avoids such 
a problem in the present approach. 
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FIG. 3. ao vs X using e(k)AF, and 5'q=o,7 = 0.018 eV. 
The inset shows experimental results for ao vs Ap, which is 
the concentration of mobile holes per CuO plane relative to 
the optimal value corresponding to YBa2Cu307 when doping 
with Pr, Ca, Zn and/or Co reduces Tc (taken from J. P. Franck 
et al., Phys. Rev. B 44, 5318 (1991); J. P. Franck et al., in 
Lattice Ejfects in High-Tc Superconductors, World Scientific 
(1992), p. 148; and references therein). 

We also remark that ao ~-{ATc/Tc){M / AM) vs x 
has the proper shape compared to experiments (Fig. 3). 
The isotope coefficient can be as low as 0.05 at Xopt- The 
minimum in ao is caused by the van Hove singularity 
in the dispersion, as remarked in previous papers. ||Tl| ] 
Away from Xopt, ao recovers the value close to 0.5 as 
in standard phononic systems. The behavior of ao is 
regulated equally by ATc and Tc i.e. ATc vs x is not 
constant in our calculation, but minimized at Xopt- 

Summarizing, here we have proposed a model where 
phononic pairing occurs between holes that are strongly 
dressed by AF fluctuations. The model may be applica- 
ble to the cuprates if the normal state correlation (^j^p 
at Xopt is proven not strong enough to produce pairing 
(currently under much discussion). Within the gap equa- 
tion formalism, we found Tc ~ 30K when the buckling 
mode of YBCO is considered. Although the hole- 
phonon coupling is much smaller than in normal metals, 
the large hole DOS boots Tc to realistic values. The 
same effect leads to an 0-isotope coefhcient that is small 
at Xopt but becomes close to 0.5 in the underdoped and 
over doped regimes. The symmetry of the condensate pro- 
duced by the buckling mode in combination with the AF 
induced hole dispersion is dj.2_y2 even in the low density 



of carriers regime. These ideas may provide a tentative 
unified explanation for several puzzling experimental fea- 
tures observed in the cuprates, especially the presence of 
an abnormal Fermi liquid at T > coexisting with a 
nonzero isotope effect. 
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